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Abstract 

The spectra of quantum dots of different geometry ( "quantum ring" , "quan- 
tum cylinder", "spherical square-well" and "parabolic" ) are studied. The 
stochastic variational method on correlated Gaussian basis functions and a 
large scale shell-model approach have been used to investigate these "artifi- 
cial" atoms and their properties in magnetic field. Accurate numerical results 
are presented for N=2-8 electron systems. 



Typeset using REVTgX 



* Corresponding author (e-mail: vargak@ornl.gov). 

^On leave of absence from Nuclear Physics Institute, Academy of Sciences of the Czech Republic, 
250 68 Rez near Prague, Czech Republic 

1 



I. INTRODUCTION 



The possibility of fabrication of artificial atoms or quantum dots with "tunable" proper- 
ties is a fascinating new development in nanotechnology. These quantum dots not only offer 
the opportunity of various applications (laser and electronic devices, memories, quantum 
gates, etc.) but they are quite intriguing physical systems in their own right. 

There are two very common ways to fabricate quantum dots. In the first method pho- 
tolithography is used to create nanoscale electrodes on the surface of heterostructures and 
the confining potential is due to the electric voltage between these electrodes Jj]. The second 
method uses material growth techniques to fabricate self- assembled quantum dots 0. 

Most of the theoretical model calculations use the effective-mass approximation to study 
the energy levels or other properties of the electrons confined in quantum dots. These 
calculations address the low-energy sector where the interband mixing is assumed to be 
negligible and the periodic crystal potential is taken into account through the effective mass 
and dielectric constant. In these models the electrons move in an external confining potential 
and interact via the Coulomb interaction. 

Given the geometry of the quantum dot and the parameters of the heterostructure, the 
confining potential can be determined by a self-consistent calculation. This is not, however, a 
trivial task and most work on quantum dots uses simple model potentials. The confinement 
is generally very strong in the vertical z direction creating quasi two-dimensional (2D) 
systems. The confinement on the xy plane is most often assumed to be parabolic. A study 
of realistic confining potentials found that this approximation is fairly good in certain cases 
but generally the confining potential might significantly differ from a harmonic-oscillator 
onei. 

The apparent similarity of "natural" atoms and quantum dots suggests the applica- 
tion of sophisticated theoretical methods used in atomic physics and quantum chemistry 
to calculate the properties of quantum dots. Parabolically confined 2D quantum dots have 
been studied by several different well-established methods: Exact diagonalization techniques 
@J5|], Hartree-Fock approximations and density functional approaches PJT0|. Quan- 

tum Monte Carlo (QMC) techniques have also been used for 2D |TT]-|T^,|T^,|T^] as well as 
three-dimensional (3D) structures [Jl4| . Few-electron artificial atoms of spherical 0] and 
cylindrical [|B| symmetry in 3D have also been investigated in variational and Hartree-Fock 
frameworks. The strongly correlated low electronic density regime got much of attention 



due to the intriguing possibility of the formation of Wigner molecules [g,^,^ ■ The novelty 
of this paper is that we use a correlated basis function which gives very accurate results for 
different few-electron systems. The accuracy becomes important when one studies subtle 
properties such as level orders, weakly bound states, etc. We investigate different models of 
quantum dots suggested by various authors. 

In this paper the variational method has been used to solve the few-electron Schrodinger- 
equation. Two different trial function sets have been applied. 

In the first case the wave function is expanded in terms of harmonic-oscillator shell model 
(SM) basis states. This basis forms a complete set and the energy is obtained by diagonaliz- 
ing the corresponding eigenvalue problem. The only approximation is the truncation of the 
basis. The dimension of the harmonic-oscillator basis quickly increases with the number of 
single-particle states included and even the powerful Lanczos method becomes unfeasible. 
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The advantage of this approach is that it is simple. Once the matrix elements are calculated 
there is no need for optimization of the basis set. In addition, we may improve on this ap- 
proach by utilizing the starting-energy independent two-body effective interaction [53| that 
takes into account two-electron correlations from the excluded space. Again, no additional 
optimization is needed as the effective interaction does not depend on any extra parameter. 

In the second approach a Gaussian basis is used. This basis is nonorthogonal and over- 
complete. The trial function depends on the parameters of the Gaussians and one has to 
optimize the parameters to get the best energy. The most adequate basis functions are 
selected by the stochastic variational method (SVM) |2"2] , |2"3| . The advantage of this basis 
lies in its flexibility. A relatively small number of basis functions give very accurate results 
provided that the parameters are carefully optimized. 

We have carefully compared the results obtained by these basis states to test the accuracy 
of the energies and other physical properties. Despite of the fact that several calculations 
exist for the 2D case, only Ref. ]TB|] reports numerical values of energy to the best of our 
knowledge. These quantum mechanical systems provide us with very good tests of different 
approaches, and therefore we think that it is important and useful to tabulate the energy 
and other quantities of the quantum dots. This may serve as a benchmark test to compare 
different methods. 

An intriguing feature of these systems is that the strength (and shape) of the confining 
potential can be changed. Unlike the natural atoms the relative importance of the pairwise 
Coulomb interaction and the external potential can be tuned. 

We have calculated the ground and the first few excited states. The ground and excited 
states are characterized by the total orbital angular momentum L and spin S. The order of 
the levels depends on the external (confining) potential. We have investigated how the level 
order changes as the parameter of the potential is varied. 

The next section introduces the basics of our formalism. The results of the calculation 
for different systems are presented in sec. III. The last section is devoted to discussion and 
summary. 



II. THE FORMALISM 

We investigate a system of N e electrons confined by the potential V con (r). The Hamilto- 
nian is 

Ne ( k 2 \ e 2 Ne 1 

*=£(-s^+w.>) +7 gi^. (1) 

In Eq. (HD, m* is the effective mass of the electron, and e is the dielectric constant of the 
semiconductor. In the following (if not explicitly specified otherwise) we will use effective 
atomic units, defined by H = e 2 /e = m* = 1. In this system of units, the length unit 
is the Bohr radius (a = h 2 /m e e 2 ) times e/(m*/m e ), and the energy unit is the Hartree 
(H = m e e 4 /h 2 ) times (m*/m e )/e 2 where m e is the mass of the electron. For the GaAs dots 
we consider here, e = 12.4 and m* = 0.067m e , and the effective Bohr radius <2q and the 
effective Hartree H* are ~ 97.94 A and ~ 11.86 meV, respectively. These effective length 
and energy units are called atomic units (a.u.) in what follows. 
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In the variational method the trial wave function is expanded in terms of basis functions: 

* = (2) 

i 

and the variational energies are obtained by solving the generalized eigenvalue problem 

- EOiJcj = 0, H id = (QilHlfy) and () tJ <I< ; ). (3) 

3 

The energy eigenvalues Ei,E 2 , ... are variational upper bounds of the energies of the ground 
and first, second,.... excited states. 



A. Correlated Gaussian basis functions 



The correlated Gaussian basis is defined in the following way: 
$<(r) = A |exp(-irAir)^ LAfL (r)x5A/ s | , 



(4) 



where A is the antisymmetrizing operator for the electrons and r = (r\, ...,r^ e ) stands for 
a set of spatial coordinates of the electrons, r^r is a short-hand notation of the quadratic 



form Ylfk=i(Ai)jkYj ' r fe) where Ai is an N e x N e symmetric positive-definite matrix whose 
elements are variational parameters. Both the spin function xsm s and the angular function 
@lMl (r) are constructed by successively coupling the corresponding single-particle functions: 



XSM S 



£i(1Xi(2) I ^(3) 

.2 2 J S12 2 



S123 



SM S 



and 



'LMj 



[y h (vi)y h (r 2 )} ll2 yi 3 (v 3 ] 



£l23 



LM L 



(5) 



(6) 



where £i m and 3^/m( r ) — r '^m(r) are the spin and angular functions of the electron. 

The Hamiltonian we consider in this paper contains no term which couples the spin and 
orbital angular momentum, and commutes with the total spin and total orbital angular 
momentum or their z components when the uniform magnetic field is applied in the z 
direction. There is no coupling between the spin and the orbital part in the basis function 
ofEq. (|). 

The correlated Gaussian function can be rewritten in a more intuitive form: 



<>X P I -9 £ a ki(*k - ri) 2 - - fori 
z k<i z k=i 



(7) 



a>ki and (3k can be expressed by the elements of A and vice versa. The advantage of this 
notation is that it explicitly connects the nonlinear parameters a^- to the pair correlation 
between the particles i and j and thus explains the name "correlated Gaussians". The 
second part, exp(— \ J2k=i Pk r 1)i 1S a product of independent single-particle Gaussians. 
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B. Separation of the relative and the center- of- mass motion 

When a system is subjected to an external field, its relative and center-of-mass motion 
cannot be separated. The harmonic-oscillator confinement is the only exception. To separate 
the relative and center-of-mass motion one can introduce a relative (e.g. Jacobi) coordinate 
system (xi, ...,Xjv e ) and rewrite the Hamiltonian and the wave function in terms of the 
relative coordinates. If the Hamiltonian is translational invariant, then the center-of-mass 
Hamiltonian H cm = -^V^ + \Mu 2 R 2 (M = Y,i=x ™>* and R = is the center-of- 
mass coordinate) can be separated. The eigenenergies of the center-of-mass Hamiltonian 
are E uX = (2u + A + 3/2)hu. In the following we are interested in the energies e int of the 
H mt — H — H cm internal Hamiltonian. The lowest energy states of the system are given by 

If the center-of-mass and the relative motion can be separated, then we use Eq. (7) by 
setting fa = and the angular part (Eq. (6)) is replaced by a similar expansion, but now 
the arguments of the spherical part are the relative coordinates Xj. In that case we have 
only N e — 1 independent variables. 

For the pure harmonic-oscillator confinement case one can use both the relative (the 
center-of-mass separated) and the single-particle coordinate (the center-of-mass motion in- 
cluded in the Hamiltonian and the wave function) approach. The energy converges to the 
same value (except for the trivial E 00 shift). The solution in the relative coordinate approach 
is, however, much easier as convergence is much faster because the center-of-mass degrees 
of freedom is decoupled. 



C. Stochastic variational method 

The energy crucially depends on the variational parameters. The optimal nonlinear 
parameters are selected by the stochastic variational method |2^j2^]. In each step of this 



procedure, JC different Ai are generated by randomly choosing the values of atki and fa from 
the [0, j3] interval. The parameter set which gives the best variational energy is selected and 
the function corresponding to that parameter set is added to the set of basis functions. The 
trial function also depends on the intermediate coupling quantum numbers (sx2, S123, ...) and 
(lx, I2, 1x2, h, ■••)■ These possibilities are also randomly tested during the optimization of the 
basis. 

Our stochastic selection procedure uses the following steps: 

(1) To set up a new basis or enlarge an existing one: 
Let us assume that the basis set has M — 1 elements. One generates /C random basis states 
and calculates the energies Ej^i (i = 1,...,/C) with the new A/"-dimensional bases which 
contain the iih random element and the preselected M — 1 basis elements. The random state 
which gives the lowest energy is selected as a new basis state and added to the basis. The 
variational principle ensures that the energy of the A/"-dimensional basis is always lower than 
that of the M— 1 dimensional one. This procedure therefore guarantees to lead to a better 
and better upper bound of the ground state energy. Notice that as the Af — 1 dimensional 
basis is orthogonalized this method does not require the diagonalization of A/"-dimensional 
matrices |22| , |23|| . The energy gain, = E^ — E^-x, shows the rate of convergence. A 



calculation of good convergence gives ~ 0. 
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(2) Refinement: To improve the energy of a basis: 

In the previous step only the newly added element is optimized, but the rest of the basis 
is kept fixed. In the refinement we keep the dimension of the basis fixed and try to replace 
the kth basis element with /C randomly generated elements. If the best energy obtained by 
substituting the kth basis state with the random candidate is lower than that of the original 
basis, then the kth basis state is discarded and the new random state is included in the 
basis. This procedure is cyclicly repeated for k — 1, . . .,jV. As the dimension of the model 
space is fixed, this step does not necessarily give lower energy, but in practice in most cases 
it does. Actually if one cannot find better basis elements, that is an indication of a well 
converged energy/basis. Again no diagonalization is needed in this step when starting from 
an orthogonalized basis. 

(3) Optimization by "fine tuning" of the parameters: 

In step (2) the parameters are randomly selected irrespective of their previous values. This 
certainly helps to avoid the traps of local minima, but if one is already (presumably) close to 
the "global" minimum then the chance to move closer to it is small. If the basis parameters 
are "reasonably" optimized or further repetition of step (1) or (2) does not lead to appreciable 
changes, one may try to change the basis parameters by selecting new parameters in the 
vicinity of the existing ones. That increases the probability of finding the nearby minimum. 
In practical calculations this step was implemented by requiring the new random parameters 
to be in the [0.8a, 1.2a] interval (a is the previously chosen parameter). In this case the 
basis optimization is done in exactly the same way as in step (2). The only difference is that 
the search interval is limited and defined by the previous parameters. 

A combination of steps (l)-(2)-(3) is repeated until the required accuracy is reached. A 
practical and economical way to set up a basis is to generate M elements (J\f = 20 or 40 
is a reasonable choice) by using step (1). Then repeat step (2) for each basis state several 
(say 3-5) times. Use step (1) once more to enlarge the basis by adding M elements to it 
and repeat step (2) as described before. After reaching a certain basis size where further 
repetition of steps (1) and (2) does not yield considerable improvement then try step (3). 

This basis selection procedure proved to be quite reliable and provides a very accurate 
solution. More details can be found in |22| . 



D. Harmonic-oscillator basis 

Alternatively, we also set up a harmonic-oscillator basis ||25||. In this case the basis 



functions are 

$ i (r) = det{(^. m .(r i )^0'))}, (8) 

where the single-particle function ip nj i jmj is a harmonic-oscillator function. This basis de- 
pends on only one parameter, the harmonic-oscillator width. For harmonic-oscillator con- 
finement this is chosen to be equal to the oscillator frequency of the potential. In this way 
the harmonic-oscillator basis functions are eigenfunctions for a noninteracting system. 

This is an orthogonal basis and the Hamiltonian matrix is sparse. The Lanczos method, 
in particular the Many-Fermion Dynamics shell-model code [p6fl , is used to find the lowest 
eigenvalues. In the diagonalization we used all states up to Y^f=\{^ n i + h) < 



' max • 
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The basic difference between the two bases is that the Gaussian basis is explicitly cor- 
related. It explicitly depends on the — r,| distances, so it is better suited to describe 
the electron-electron correlations. At the same time the harmonic-oscillator basis is simpler 
because no optimization is needed. 

An advantage of the harmonic-oscillator basis is the fact that we may alternatively per- 
form the calculations in the Jacobi coordinates with the center-of-mass degrees of freedom 
removed. It is straightforward, although numerically intensive, to construct an antisym- 
metrized harmonic-oscillator basis depending on the Jacobi coordinates ||24|| . Depending on 
the problem, we may choose the more efficient basis. For N e = 3,4,5 electron systems it 
turns out that the use of Jacobi coordinates is more profitable. For larger number of elec- 
trons, it is more efficient to use the single-particle coordinates and the Slater determinant 
basis (||). 

As the harmonic-oscillator frequency is fixed as described above, the only parameter of 
the calculation is the model space size characterized by A max . In the present calculations 
we use as large N max as possible, typically iV max = 15 — 33 for N e < 5 and iV max = 8 — 12 
for larger systems. 

A speed up of convergence can be achieved by utilizing the effective interaction approach 
that was succesfully applied in the ab initio shell model calculations for few-nucleon systems 
and light nuclei [ 2"! . While it is crucial in the nuclear physics application to use the 



effective interactions, in the present electron systems the effective interaction provides only 
minor improvement. In some cases, however, it brings the SVM and the SM results to much 
closer agreement. The details of how the effective interaction is computed from the bare 
Hamiltonian, here the harmonic-oscillator and Coulomb interaction, is given, e.g., in Refs. 



24j ) |25| . The basic goal of the effective interaction is to take into account, in this case two- 
electron, correlations from the excluded space, i.e., from the space containing excitations 
above A max . A formal difference from the nuclear case is that here the harmonic-oscillator 
potential is a real binding potential, while in the nuclear application it is a model potential 
representing the mean field formed by all nucleons, which is added/subtracted to the real 
nucleon-nucleon interaction in order to facilitate the effective interaction calculation. 

We note that when the m-scheme basis @ is used the good quantum numbers are checked 
by evaluating the mean values of relevant operators, e.g., J 2 , L, and S for each eigenstate. 



E. Magnetic field 

In external magnetic field the kinetic energy operator is replaced by 

1 :P 2 - A(P* + -A,) 2 . (9) 



2m* * 2m* c 

We consider a uniform magnetic field B = (0, 0, —B). By taking Aj = — o r i x B the above 
expression can be rewritten in a more detailed form 

-^(P* + - c ^f = ~^ 2A * + \m*{u c /2f{x^ + yj) - \u c l zi , (10) 

where l Z i is the z component of the orbital angular momentum of the ith electron. The 
cyclotron frequency for the parameters we use in this paper reads as 



7 



f k o c = — = ^p B B = 0.145725 (iT), (11) 

m*c m* 

where the Bohr magneton is p B = eh/(2m e c) = 0.05788 meV/T. The interaction of the 
magnetic field with the spins leads to the Zeeman term, —g* p B Bs Z i, where s Z i is the z 
component of the spin of the zth electron and g* is the effective ^-factor of the electron. 
The Zeeman term leads to the splitting of the energies for different spin orientations. As 
the Hamiltonian with this term still commutes with the z component of the total spin, 
S z = Si=i s Z i, the energy shift is simply given by —g*p B BS z and one can easily add this to 
the energies presented in the following. This energy is not included in what follows. 

The correlated Gaussians defined above are not ideally suited for systems in magnetic 
field, because the basis functions belonging to different orbital angular momenta would be 
coupled by the Hamiltonian. This coupling would require an infinite series of orbital angular 
momentum states, which is obviously out of question. To avoid this, we choose a deformed 
form of the correlated Gaussians (DCG) |27f : 



i I N e ^ Ne | 

ex P \ ~2 ^ Ai i pi 'Pi~ 2^ BijZiZj I ' ^ 

where the nonlinear parameters are different (and independent) in the xy and z directions 
(pi = (xi,yi)). This extension brings a great deal of flexibility by allowing a separate 
description on the xy plane and along the z axis. The Hamiltonian does not commute with 
L 2 but it does with L z . The eigenfunctions have good quantum number M of L z . Note 
that we will use M for the orbital angular momentum quantum number in 2D and L for the 
one in 3D. The above form of the DCG belongs to M — 0. To allow for M^fl states we 



multiply the basis by [22 



Hum, (is) 



where 



t ( \ _ / ( x + ^y) m l° r nonnegative integer m , > 

\ (x — iy)~ m for negative integer m. 

Thus our variational basis function reads as 

$ M (r) = A I WPi)^ ex P ^ - \ Yl A vPi ■ P3~\Yu B U~-~-.?j | > (I 5 ) 

where M = mi + rri2 + ... + mjv e . 

The above basis is defined for 3D cases. It is used not only for solutions in the presence of 
magnetic field but also for external potentials with cylindrical symmetry. For 2D calculations 
the same form is used except that the third component of the vectors are dropped (or 
equivalently By = is assumed). 
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III. CALCULATION 



A. Harmonic-oscillator confinement in 2D 



The harmonically confined 2D systems received much theoretical attention and this is a 
very good test case to gauge the accuracy of different approaches. In this case the confining 
interaction takes the simple V coa (r) = ^m*u 2 r 2 form. The single-particle energy of the 
harmonic-oscillator potential is given by (2n + |m| + l)hu, where n = 0, 1, 2, . . ., and m = 
0, ±1,±2, .... In Table I we compare our results to the "exact diagonalization" and 
the QMC methods [] I - |T3|j for the N e = 3 electron system. We have carefully optimized the 



parameters and repeated the calculation several times to check the convergence. Our result is 
expected to be accurate up to the digits shown in Table I. In principle the QMC calculations, 
except for the statistical error, give the exact energy of the system. In practical cases the 
famous "minus-sign problem" forces the QMC approaches to use certain approximations 
(in Ref. [13,12] the "fixed-node" method has been used). The slight difference between our 



results and the QMC values is probably due to this fact. The energies for both the ground 
and excited states are in good agreement. Our results are slightly better than the other 
calculations in each case. 

In Table II a similar comparison is presented for N e = 2 — 6 electron systems. The QMC 
results |TB| quoted in Table II are obtained by very careful calculations and their statistical 



error is very small. Note that the confining strength is slightly different in the calculations 
presented in Tables I and II. This table also includes the virial factor 

N e 

r ] = 2(T)/(W), (W) = (£r i -V i V ist ), (16) 

i=i 

where V\ nt is the "interaction part" of the Hamiltonian, including the confining and the 
electron-electron interactions. The virial factor is unity for the exact wave function. 

Our result is in excellent agreement with the QMC predictions |13|] in all but one case 



(iV e = 4). The QMC renders the (M, S) = (0,0) state to be the ground state and the 
(M, S) = (0, 1) state to be the first excited state, which is a violation of the Hund's rule. 
The shell filling and Hund's rule has been experimentally investigated in Ref. PS and it is 



found that a circular dot obeys the Hund's rule. According to the Hund's rule the ground 
state of a system with a well developed shell structure is in the maximum spin state allowed 
by the Pauli principle. See the Appendix for an example of the N e = 4 electron case. The 
violation of the Hund's rule in this system is also observed in another QMC calculation ||11|| . 
This latter calculation predicts a relatively large energy difference between the (0,0) and 
(0,1) levels, but it is somewhat less accurate for N e = 4, using only the lowest Landau levels. 

Our calculations, in agreement with the Hund's rule, predicts the (M, S) = (0, 1) state 
to be the ground state and the (2,0) and (0,0) states to be the lower excited states. This 



contradicts the results of [P] , [l3|] but is in agreement with the other QMC calculation [12 
Our other energies are very close to the QMC results: The agreements for N e = 5 and 
N e = 6 electron systems are very impressive. 



We define the pair correlation function 

2 

N e (N e - 1 



p ( r > r o) = aFtw tt <*I E 5 ( r * - R - r )*fo - R - r o) l*>- ( 17 ) 



l<3 
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Here r is a fixed vector and its magnitude is chosen to be equal to (^\ J2i l r « — R-l \*&)/N e . The 
function P(r, tq) gives us information on where one electron located at ro sees other electrons. 
Figures 1 and 2 display the pair correlation functions for the ground state (M, S) = (1, 1/2) 
and the first excited state (M,S) = (2,3/2) of N e = 5 electron system. Both figures 
show qualitatively similar features. For u — 1, the confinement potential is strong and 
the contribution of the single-particle energies to the total energy is larger than that of 
the Coulomb potential. The electrons are confined in a rather compact region so that the 
contour map does not show clear four peaks. On the contrary, for u — 0.1 the effect of 
the confinement becomes weak and the contribution of the Coulomb potential is larger than 
that of the harmonic-oscillator part. The size of the system grows and we see clearly well- 
separated pentagon-like structure. 

Next we present in Table III an example where the magnetic field is nonzero. Again, the 
energies are in good agreement with the QMC [11] and diagonalization || methods. In 2D 



the inclusion of the magnetic field leads to a change of the harmonic-oscillator frequency 



CO 



^ + (u c /2) 2 (18) 



and an energy shift by — ~Mhuj c , so we expect that our results are as accurate as those for 
the zero field case. The accuracy is also indicated by the virial factor included in Table III. 

We have improved the prediction of the diagonalization method M. The diagonalization 
method would give the "exact" solution in infinite model space. In practice the diagonaliza- 
tion is always limited to finite dimensions. The slight disagreement between our and QMC 
results might be due to the statistical (and/or fixed node) error of the QMC calculation. 

The states listed in Tables I and III (i.e., (1,1/2), (2,1/2), (3,3/2)) are quoted as the 
lowest-energy states in several papers P JIT|J12]| . Figure 3 shows that this is not the case. 
The level order at B = is (1, 1/2), (0, 3/2), (2, 1/2), (0, 1/2), (3, 3/2). Figure 3 shows that 
the spin unpolarized (1, 1/2) state is the ground state in the weak magnetic field limit and the 
(2, 1/2) unpolarized state becomes the ground state in a very small interval of the magnetic 
field strength. The spin polarized states become the ground state above B =2.5 T. The 
figure also shows that the lowest spin polarized state is the (0,3/2) state for weak field. 
For stronger field the (3, 3/2) and then the (6, 3/2) states become the lowest spin polarized 
(and ground) state, following the (3, 6, 3n) "magic" sequence. Other spin polarized states 
(e.g., (1,3/2) etc.) never become the lowest state. The explanation of the magic sequence is 
very simple. In the spin polarized case all electrons have to occupy different orbits. As the 
magnetic field gets stronger, the single-particle states belonging to positive orbital angular 
momentum quantum numbers (mj = 0, 1, 2, 3, ...) are energetically more favorable than those 
with negative ones. The M = 3 state ((m 1 ,m 2 ,m 3 ) = (0, 1,2)) is therefore lower than the 
M = 2 state (which requires (0, —1,3) or (1,2, —1), etc.). For the weak magnetic field the 
above argument does not hold in general and the lowest polarized state is M = with the 
(0, 1, -1) orbits. 

Similar picture is valid for N e = 4 (see Fig. 4). In the very weak field regime the 
unpolarized (M, S) = (0, 1) state is the ground state. By increasing the magnetic field, 
the spin polarized M = 2 state ((mi, m 2 , m^, m^) = (0, 1, —1, 2)) becomes the ground state 
before the "magic" M = 6 (0, 1, 2, 3) state takes over. 

Figure 4 reassures that the (M, 5") = (0, 1) state is the ground state and the (M, S) = 
(0, 0) is an excited state for zero magnetic field: Both states belong to M = 0, and therefore 
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the change of the magnetic field simply changes the harmonic-oscillator frequency (see. Eq. 
(18)). The figure thus shows that the order of these two states remains the same for different 
harmonic-oscillator frequencies. 



B. Harmonic-oscillator confinement in 3D 



We have calculated the energies of the ground and first few excited states of 3D few- 
electron systems confined by a harmonic-oscillator potential (V con (r) = ^m*uj 2 r 2 ). The 
results for different values of the oscillator frequency are compared in Tables IV-X. All 
intermediate spin coupling possibilities (si 2 , S123, •••) are included in the trial function. The 
partial wave components (Zi, l 2 , ■■■) are included up to J2i=i X k < 6. 

The two-electron case is relatively simple and it is analytically solvable for certain fre- 
quencies [29]. For uo = 0.5 for example, the exact energy is 2 a.u. [29] and we can easily 
reproduce this value up to several digits as shown in Table IV, where the energies of other 
low-lying states are also listed. Three very different oscillator frequencies are used to test 
the accuracy of the method under different circumstances. In the case of u = 0.01 the 
confinement is extremely weak and the Coulomb interaction governs the dynamics. In the 
other limiting case the confinement is very strong (a; = 10). Another reason for choosing 
these values is that we want to study the ordering of the energy levels as a function of the 
strength of the confining interaction. In the two-electron case, for example, there is a level 
crossing between the state (L, S, n) = (3, 1, — ) and the first excited state of (0, 0, +). The 
order of the other states listed in Table IV does not change. 

The energies of the ground and excited states calculated by the correlated Gaussian and 
the harmonic-oscillator shell model basis are compared for N e = 3 — 6 electron systems in 
Tables V-VIII. Both methods give very similar results: The agreement is especially good 
for the ground and first excited states. For higher excited states the Gaussian basis gives 
slightly less accurate energies because it is significantly more difficult to optimize the basis 
for excited states. In addition to the present results we note that the harmonic-oscillator 
calculation dependence on the model space size for N e = 3, 4 and u = 0.5 was discussed and 
tabulated in Ref. []30| . 

The agreement is especially striking for N e = 3. Almost all digits are equal for most 
of the calculated cases. It is interesting to compare the order of the states in 2D and 3D. 
In the 2D case for N e = 3 the energy levels of the first few states follow the order of 
(M, S) = (1, 1/2), (0, 3/2), (2, 1/2), (0, 1/2), while in 3D the levels are ordered as (L, S, tt) = 
(1, 1/2, -), (1, 3/2, +), (2, 1/2, +), (0, 1/2, +). This shows that the lowest levels are built up 
from the same single-particle states: In the 3D ground state two electrons are in the I = 
orbital and one is in the I = 1 orbital. The first excited state has two electrons in the 
I = 1 orbital, which are coupled to L = 1 because their spin must be parallel to build up 
S = 3/2 with the third. In the 2D case they are in the m = 1 and m = — 1 orbitals and 
their total orbital angular momentum is M = 0. The higher excited states have similar 
correspondence. The same similarity occurs for N e = 4 ((M, S) = (0, 1), (2, 0), (0, 0) in 2D 
and (L, S, n) = (1, 1, +), (2, 0, +), (0, 0, +) in 3D). For example, the 3D ground state has two 
electrons in the I = and two electrons (with parallel spin) in the I = 1 orbital and the 
two electrons in the I = 1 orbital are again coupled to L — 1. In the 2D ground state the 
two electrons are in the m — 1 and m — — 1 orbitals and the orbital angular momentum 
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is M = 0. With respect to the single-particle state occupations there is of course a big 
difference between the 2D and 3D case. In 2D the shell fillings occur at N e = 2, 6, 12, 20, . . . 
etc., while in 3D the shells are filled at N e = 2,8,20,40, .... For N e > 6 particle systems 
the single-particle components of the wave functions in 2D and 3D might be quite different. 

The addition energy is conveniently used to show the shell closure which occurs at a 
specific electron number. The addition energy Afx(N e ) is defined by 



where the chemical potential //(iV e ) is the increase of the ground state energy by adding one 
electron to the ground state of N e — 1 system: 



The shell or half shell closure is reflected by a sudden increase of Afx(N e ) at a certain N e or 
the change of the differential capacitance given by e 2 /Afi(N e ). This is because the electron 
needs much energy when it fills an orbit across the degenerate orbits of a shell or goes beyond 
the half-shell due to the Hund's rule. The addition energies of the harmonically confined 
electrons in 2D and 3D are compared in Figs. 5 and 6. In 2D the addition energy shows a 
large peak at N e = 2 and a smaller peak at N e = 4. The former corresponds to the filling of 
the n = 0, m = orbit, while the latter is a reflection of the half shell filling of the degenerate 
orbits n = 0,m = ±1, which can be understood by the Hund's rule. By decreasing uj the 
level spacing of the single-particle orbits becomes smaller and the correlation due to the 
Coulomb interaction takes over the shell structure. This explains why the peak at N e = 2 
disappears for u — 0.1. The behavior of the 3D addition energy is similar to the 2D case. 
One difference is that the half shell filling occurs at N e = 5 because the relevant orbit is 
I = 1 and can accommodate six electrons. 

The results for N e = 5 and N e = 6 are somewhat less accurate and the agreement 
between the SVM and the shell model is not as good as for N e < 5. The SVM seems to be 
more accurate than the shell model for week confinement, where the role of the Coulomb 
interaction is more pronounced and it is more difficult to take the Coulomb correlation 
into account with the shell model basis. At the same time it is more easy to use the shell 
model approach for larger systems (see Tables IX and X), while the SVM becomes very time 
consuming beyond N e = 6. 

A general feature of the results is that the excited states change their level orders as 
the harmonic-oscillator strength changes, but the ground state always remains the same. 
We have very carefully tested this property and we do not find any level crossings with the 
ground state. 

Other insights of the relation between the 2D and 3D systems can be gained by comparing 
the expectation values of the kinetic, confining and Coulomb operators. Tables XI and XII 
show the contribution of the Coulomb, kinetic and confinement parts of the Hamiltonian 
to the total energy. The contributions are nearly equal in the uj = 0.5 case. Just as one 
expects in the strong confinement case (uj = 10) the kinetic and confinement energies are 
strongly enhanced, and the Coulomb energy is relatively small but not negligable. On the 
other hand, in the weak confining case the Coulomb interaction dominates. 

The contribution of the confining interaction, and thus the kinetic energy is of course 
larger in 3D. If the electrons would not interact then both the kinetic and the harmonic 



A/i(i\Q = fj,(N e + 1) - fi{N e ) 



(19) 



fi(N e ) = E(N e ) - E(N e - 1). 



(20) 
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confinement energies would be 1.5 times increased in 3D compared to the 2D case. In 
the interacting case the kinetic and confinement energy increase is roughly 1.5 for N e = 2 
and N e = 3. For larger systems the increase is smaller. On the other hand, the Coulomb 
correlation energy is smaller in 3D than in 2D because there is larger space available in 3D 
for the electrons. 



C. Spherical square well 

As an alternative to the harmonic confinement one can consider a spherical square well 
model of the 3D quantum dots. In this case the electrons are confined by a square well 
potential: 

The square well potential is analytically solvable for one particle case. The eigenenergies 



E can be determined from the transcendental equation, y V — |i?|cot(y 2(V — \E\)R) = 

— yJ\E\ (for / = 0, and atomic unit used). Our SVM numerical approach virtually exactly 
reproduces the analytically determined energies. 



Spherical quantum well-like quantum dots have been studied in Ref. |fL7 |. Unlike the 
harmonic-oscillator potential, the spherical well can only hold a certain number of electrons. 
The number of electrons that a spherical well can bind depends on VqR 2 . Figure 7 shows 
the energies of few-electron systems confined by a spherical square-well potential in 3D as 
a function of the radius R. A spherical well can only bind an electron if 7r 2 /8 < V$R 2 . In 
our example Vq = 10 and therefore the one electron bound state appears when 0.35 < R. 
By increasing the radius the two, three, ... etc. electron systems may become bound in the 



well (see Fig. 7). This potential parameter has been used in Ref. [17| to simulate quantum 



dots in GaAs/Ali_ x GaaAs with x ~ 0.1. Note that due to the difference in units the radius 
used here corresponds to v^2 times the one used in Ref. fl7 ]. 



A comment is in order concerning the energy curves in Fig. 7 (and in Figs. 8 and 13 
in later subsections). If the N e electron system has a bound ground ground state then our 
calculation converges to the energy of that state. If there is no bound state in a given 
potential then the energy convergies to the lowest relevant threshold, which is in this case 
the energy of the N e — 1 electron system. In the figures the system is bound if the energy of 
the N e electron system is below the corresponding N e — 1 electron system. Strictly speaking, 
for unbound (N e electron) states the energy of the N e electron and N e — 1 electron systems 
should be equal. The convergence of the energy of the unbound N e electron system to 
the energy of the N e — 1 electron system is rather slow, so one needs many basis states to 
describe the "free" electron. Therefore the fact that the energy curves of the unbound states 
are above the corresponding thresholds is the consequence of our spatially limited basis. By 
using more basis states and by allowing them to go far outside of the range of the confining 
interaction one would get the same energy for the N e — 1 and the unbound N e electron 
systems. 

We have found no "phase transition" in N e = 2 and N e = 3 electron systems. The 



authors of Ref. flTfl have investigated the energy of the lowest spin polarized and spin 
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unpolarized N e = 2 and N e = 3 electron systems as a function of the radius of the square 
well. They found that beyond a certain radius the spin polarized state becomes lower than 
the spin unpolarized ground state. We have very carefully investigated these systems and 
have not observed this "para to ferromagnetic phase transition". The same authors in a 
later paper |18| investigated a harmonically confined two-electron system and found that 
the spin unpolarized to spin polarized transition is most likely an artifact of the neglection 
of part of the electron-electron correlation in Hartree-Fock calculations. 

We have increased the radius gradually (see Fig. 7). As the Coulomb repulsion decreases 
the energy of the system gets smaller and smaller, converging toward the energy of the 
noninteracting electrons in the quantum well. If there is no Coulomb interaction then the 
energy of the spin polarized and unpolarized electrons is the same, so by increasing the 
radius both converge to the same energy. In our present example (Vo = 10) the energies of 
the lowest lying spin polarized and unpolarized states are nearly degenerate beyond R = 15, 
but we observe no level crossing between them. 



D. Cylindrical well: "Quantum cylinder" 

In this section we present a calculation for cylindrical quantum dot. Similar case has 



been considered in Ref. [ff9| in an unrestricted Hartree-Fock framework. The confinement is 
defined as 

Kon(r) = (-K> if (x' + y 2 ) 1 / 2 <Rand\z\<a (22) 
1 otherwise. 

In this case the spherical symmetry is broken, and only the z component of the orbital angu- 
lar momentum is conserved. We have to use the DCG basis functions that were introduced 
for magnetic field. 

First we consider a model potential with Vo = 10 and R = 1 and change the "thickness" 
of the dot (the height of the cylinder) from a — 10 to a — (in a.u.). In this way we 
transform the system from a rod-like (a = 10) geometry to a 2D disk (a = 0). Just like in 
the case of the spherical quantum well a quantum cylinder can bind only certain number of 
electrons, depending on the potential parameters Vo, R and a. The energy dependence on 
the thickness of the cylindrical dot for N e = 1 — 4 electron systems is presented in Fig. 8. 
The figure shows that as one expects the cylinder can hold more and more electron as the 
size (in our case the height) increases. The real interesting thing here is that the order of 
the energy levels also depends on the height of the cylinder. For long, rod-like cylinders the 
ground state tends to be the M = orbital angular momentum state. This probably means 
that the electrons are equidistantly positioned along the z axis. By decreasing the height we 
approach a disk-like geometry which is somewhat similar to the 2D harmonic confinement 
discussed earlier. And indeed, the level order changes (see Fig. 8) and one has the same 
level order as in the 2D harmonic confinement case. In this way we have found an interesting 
transition: By decreasing the height of the cylinder the (0, 1/2) and (1, 1/2) N e — 3 electron 
(also (0,0) and (0, 1) N e = 4 electron) ground/excited states change their order. 

Figures 9 and 10 show the density of the N e = 2 electron system in a "long" cylinder 
(a = 10) along the radius and along the symmetry axis, respectively. The radial density 
distribution of the triplet (1, 1) and singlet (0, 0) states are very similar. Both peaks around 
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0.5 a.u. and the tail goes a little bit outside of the cylinder. The density of the triplet and 
the singlet state in the z direction, however, are very different as shown in Fig. 10. The 
cylinder is so long that the two electrons can be far away from each other to minimize the 
Coulomb repulsion and the density tail hardly goes outside the cylinder. 

Figures 11 and 12 show the density distribution for a "thin" cylinder (a = 1). Decreasing 
the height of the cylinder the triplet state becomes less and less bound. The radial density 
distribution of this very weakly bound triplet state has one peak in the cylinder and another 
outside. The density distribution of the triplet pair along the z axis has a peak inside the 
box but also shows a long tail outside. So one of the electrons is inside the box and the 
other one is already mostly outside. By decreasing the height further the electron which is 
outside will be unbound, and the cylinder will be able to bind only one electron. The singlet 
state is still bound for a = 1 but it will also be unbound if one decreases the height of the 
cylinder further. 



E. "Quantum ring" 



Ring-like nanostructures have been grown by electron-beam lithography |3l |. The elec 



tronic and magnetic properties of a single electron quantum ring have been studied in Ref. 



We restrict our attention to a pure 2D case. An additional confining interaction in the 
z direction would cause no extra difficulty in our approach. The confinement in this case is 
defined as 



Kon(p) = { -V n < p < r 2 (23) 




This describes a square well potential in a ring between r\ and r 2 on the xy plane. 

The number of electrons bound in a ring-like potential, similarly to the previous cases, 
depends on the parameters (Vo,ri,r2) of the potential. An example of the energy levels 
in the model potential is presented in Fig. 13. The maximum number of electrons this 
potential can bind is N e = 4. In the four electron case the lowest state is the (M, S) = (0, 1) 
state just like in the case of 2D harmonic confinement. In the present model potential the 
first excited state (0, 0) is not bound with respect to the three-electron threshold. 

The density distribution of an electron in a quantum ring is shown in Fig. 14. The 
electron is along the ring between r\ and r 2 . For the potential strength used in this example 
the distribution is "wider" than the width of the ring. By increasing the potential strength 
the "inner" and "outer" tails of the density are pulled more and more inside the quantum 
well. 

This density distribution can be easily manipulated by a perpendicular magnetic field. 
The magnetic field acts as a confining harmonic-oscillator potential on the xy plane. By 
increasing the strenght of the magnetic field, the density distribution starts to move inward 
as shown in Fig. 15. In a certain very narrow region of the magnetic field strength it has two 
peaks: An outer peak centered in the ring and an inner peak which is inside the harmonic 
confinement induced by the magnetic field. If the magnetic field is stronger than a given 
value then the electron moves inside of the harmonic confinement. This geometry gives us a 
possibility of moving the electron from one well-defined position to another one by switching 
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on and off the magnetic field. Notice that we have two sharply separated peak positions in 
this case. In the case of the previously studied harmonic or spherical square well confinement 
the density distribution also moves toward the origin in the magnetic field. In that case, 
however, what we see is more like a "shrinking" of the density distribution on the xy plane. 
The change is more drastic in the case of the quantum ring. The peak of the distribution 
shifts from one position to another. 



IV. DISCUSSION 

The properties of artificial atoms created by confining electrons in quantum dots of 
different geometry are qualitatively very similar. The electrons occupy the single-particle 
orbits defined by the confining interaction. The occupancy is determined by the Pauli 
principle and the minimization of the Coulomb energy. In different confining potentials the 
energy levels are different but the basic features are very similar. 

The confining interactions considered in this paper depend on one or more parameters 
(harmonic-oscillator width, radius and strength of square- well, etc.). We have studied the 
dependence of the energy levels on these parameters. One intriguing property that we have 
found is that in spherically symmetric systems the ground state remains the same for any 
values of the parameters. Its energy level does not cross with that of the excited states. At 
the same time, the order of energy levels of the excited states frequently changes depending 
on the parameters of the confining interactions. The change of symmetry of the ground 
state of a spherical quantum well has been reported in Ref. |17|]. We have investigated 
few-electron systems in spherical quantum wells of different parameters but we have not 
observed any similar change. This confirms that such change of energy levels might be an 
artifact of Hartree-Fock Calcualtion [fTB |. Our calculation predicts that the ground state is 



in accordance with the Hund's rule for any parameter values and there is no transition from 
spin unpolarized to spin polarized states. 

The change of the ground state would give us an interesting possibility: In a two electron 
system, for example, the ground state is a spin singlet, and the first excited state is a spin 
triplet. This two-state system may serve as a "qubit", an elementary gate for a quantum dot 
quantum computer. One would prepare a dot with singlet ground and triplet excited state 
and a second one, with different geometry, where it is the other way around. The electrons 
can be moved from one dot to the other by an external electric field, for example, switching 
from S — to S — 1. The calculations show, however, that no matter how do we change 
the geometry, the ground state does not change for spherically symmetric systems. 

If the spherical symmetry is broken, for example by a magnetic field or by a cylindrically 
symmetric confining potential, then the ground state and excited state energy levels may 
cross each other. The fact that the magnetic field changes the order of energy levels has been 
studied in many works. In this paper we have presented a method which very accurately 
predicts the level crossing as a function of the strength of the magnetic field. 

The cylindrical quantum dot shows a very interesting example where the order of energy 
levels depends on the height of the quantum dot. In a disk-like cylinder the ground state 
obeys the Hund's rule, but for a longer cylinder, typically when the height of the cylinder 
is larger than its diameter, the Hund's rule is violated. It would be interesting to look for 
experimental evidence showing that the ground state of a cylindrical three- (four) electron 
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quantum dot is a (M,S) = (1,1/2) ((M, S) = (0,1)) state if the height of the cylinder is 
small and the ground state becomes a (M, S) = (0, 1/2) ((M, S) = (0, 0)) state by increasing 
the height of the cylinder as predicted here. 

Harmonically confined electron systems in 2D have attracted enormous attention. In this 
work we have also calculated 3D electron systems in harmonic confinement. The comparison 
of the 2D and 3D cases shows the effects of the quantum-well confinement in the z direction 
in quantum dots. For the same harmonic-oscillator strength the electrons are somewhat 
farther from each other in 3D than in 2D, resulting in a smaller Coulomb energy in 3D. The 
energy difference between the 2D and 3D geometry is dominantly due to the confinement and 
the kinetic energy. The qualitative features of the 2D and 3D systems are very similar in the 
case of the few-electron systems investigated here. One can make an easy correspondence 
between the orbital and spin quantum numbers of the energy levels in 2D and 3D. The 
applicability of our method is tested by calculations for very different confining strength. 
The accuracy is slightly worse for the weak confining region where one needs more basis 
functions to achieve convergence. But the accuracy is fairly good as one can judge it by the 
virial factor and by comparing with the results of different methods. In the weak confining 
region (u = 0.01, see Tables XI and XII) the contribution of the kinetic energy is fairly small 
compared to that of the Coulomb and confining interactions. This suggests the existence of 



Wigner crystal like structure both in 2D and 3D. Contrary to the prediction of Refs. ||TT]JT3 
we find that the ground state of N e = 4 system obeys the Hund's rule. 

We have also investigated an example of a ring-like quantum dot in magnetic field. This 
geometry offers an interesting possibility. In the case of zero magnetic field the electrons are 
distributed along the ring. By applying the magnetic field perpendicularly to the plane of 
the ring the electrons can be completely moved from the ring to the vicinity of the origin. 
Thus one may have the electrons in two very well separated regions. 

The major difference between the harmonic and the square well confinements (cylindrical, 
spherical and ring) is that the harmonic case can bind any number of electrons. The number 
of electrons bound in the square-well case is finite and strongly depends on the parameters 
of the potential. In that case one can predict how many electrons can be bound in a certain 
quantum dot, that is, the "capacity" |33[ of the dot. This is expected to be a more realistic 
model of quantum dots. 

The concrete potential parameters and the potential itself can only be determined ex- 
perimentally. In this work we have tried to follow the "experimentally inspired" and widely 
used potentials and parameters. The aim of this work was to demonstrate the wide range 
of applicability of the method and the investigation of different properties of various arti- 
ficial atoms. The direct comparison to experiments may not be suitable in this level. The 
general features of experimental findings may be reasonably well described by the potential 
models considered here. In the present level one assumes a Hamiltonian which models the 
quantum dot and we try to solve this well-defined quantum mechanical problem in a careful 
manner. There are of course many things which may limit the applicability of our model 
Hamiltonian, but some of the phenomena that are experimentally observed can be under- 
stood by such model calculations and hopefully some of the predictions of such models can 
be experimentally observed. 

The accuracy presented here is very useful and important in the weakly confined (but 
strongly correlated) regime where otherwise it is difficult to predict the ground state, etc. 
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One should also mention that the comparison of various methods for these quantum me- 
chanical problems greatly helps the test and developement of different quantum mechanical 



many-body approaches. A nice example can be found in Ref. [30], where the solution of 
few-electron quantum dot problems helps to test the Faddeev method which was developed 
for nuclear few-body systems. 

In summary, we have presented a large scale variational approach to describe the spectra 
and other properties of artificial atoms. Different (parabolic, cylindrical, spherical and 
ring-like) confining interactions have been investigated. The effects of magnetic field have 
also been studied. One of the aim of this paper was to introduce the method and tests 
its capabilities on various models of quantum dots used in the literature. Future work to 
investigate double quantum dots is under way. 
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APPENDIX 

The aim of this appendix is to understand in a simple 2D model how Hund's rule comes 
about for four-electron system. Let us assume that the shell spacing is fairly large, so that 
we may put two electrons in the lowest n = 0, m = orbit and other two electrons in the 
next n — 0, m — 1 and n = 0, m = —1 orbits. A properly antisymmetrized wave function is 

* = 4^W-i( 2 ) ± 4>-i(l)MmSM s ), (24) 



1 11 

tpi = — ^exp(--p 2 )(x + iy), = -^exp(--p 2 )(x - iy). (25) 

'71 I \/TC I 



where 



The plus sign in Eq. ( p4[ ) is for S = (singlet) and the minus sign is for S = 1 (triplet). These 
two states have the orbital angular momentum M = but differ in their spin configurations. 
The Coulomb energies of the two states are 

(*h " ;\^) = D±E, (26) 

\Pl- P2\ 

where D is the direct interaction energy and E is the exchange interaction energy. Whether 
the singlet state is lower than the triplet state or not depends on the sign of E. They are 
given by 

D = (^(1)^(2)1^^1^(1)^(2)) = H7| 
I Pi - P2\ 16 V 2 

E = (^(1)^(2)1^—^1(1)^1(2)) = > 0. (27) 

\Pi - Pi\ 16 V 2 
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Since E is positive, the triplet state is lower than the singlet state. 

When two electrons are put in n = 0, m = 1 orbit, they have the orbital angular mo- 
mentum M=2 and must be in the spin singlet state because their spatial part is symmetric. 
The Coulomb energy of this state is given by 

(tflp-L-Jtf) = (^(1)^(2)1^—1^(1)^(2)) = D, (28) 

\Pl — P2\ \Pl — P2\ 

which is equal to the average of the energies of the singlet and triplet states with M — 0. 

In the above discussion the kinetic energy and the harmonic-oscillator confinement energy 
are not considered, but their contributions are the same in the three states. Also the Coulomb 
interactions between the m = and m = ±1 orbits give the same contribution in these three 
states. Therefore we can conclude, in accordance with the Hund's rule, that the triplet state 
of (M, S) = (0, 1) is lower than the singlet state of (M, S) = (0, 0) and also the state with 
(2, 0) is in the middle between them. This order is exactly the same as the one we obtained 
by the more realistic calculation (see Table II). 
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FIGURES 




FIG. 1. Pair correlation function of the ground state (M,S) = (1,1/2) of 2D five-electron 
system as a function of the frequency oo of the harmonically confining potential. The white cross 
denotes i-q. Atomic units are used. 
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FIG. 2. Pair correlation fucntion of the excited state (M,S) = (2,3/2) of 2D five-electron 
system as a function of the frequency uj of the harmonically confining potential. The white cross 
denotes tq. Atomic units are used. 
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FIG. 3. Energies of the harmonically confined (Two = 3.37 meV) lowest spin unpolarized 
(S=l/2, thin solid line) and spin polarized (5=3/2 thick solid line) three-electron states in mag- 
netic field. The orbital angular momentum M of the state is indicated by the number next to the 
curve. The Zeeman energy is not included. 



22 




B(T) 



FIG. 4. Energies of the harmonically confined (Tito = 3.37 meV) lowest spin polarized (5 =2, 
thick solid line) four-electron states in magnetic field. The orbital angular momentum M of the 
state is indicated by the number next to the curve. The two thin solid curves are the 5=0 and 
5=1 (these 5 values are indicated next to the thin curves) states belonging to M=0. The Zeeman 
energy is not included. 
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FIG. 5. Addition energy of harmonically confined electrons in 2D as a function of the electron 
number, to is the frequency of the confining potential. 
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FIG. 6. Addition energy of harmonically confined electrons in 3D as a function of the electron 
number, to is the frequency of the confining potential. 
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FIG. 7. Energies of N e = 1 (solid line), N e = 2 (dotted line) and N e = 3 (dashed line) electron 
systems in a spherical quantum well as a function of the radius of the well. Lower dotted line: the 
ground state (0,0,+); upper dotted line: the excited state (1,1,—); lower dashed line: the ground 
state (1,1/2,—); upper dashed line: the first excited state (1,3/2,+). 
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FIG. 8. Energy of N e =l (solid line), N e =2 (dotted line), 7V e =3 (dashed line) and N e =A (long 
dashed line) electron systems in a cylindrical quantum well as a function of the height of the 
cylinder. (Vq = 10 and atomic units are used.) 
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3.0 




r(a.u.) 

FIG. 9. Singlet (solid line) and triplet (dotted line) radial density distribution of two electrons 
in a cylindrical quantum dot (a = 10). 
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FIG. 10. Singlet (solid line) and triplet (dotted line) density distribution along the z direction 
of two electrons in a cylindrical quantum dot (a = 10). 
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FIG. 11. Singlet (solid line) and triplet (dotted line) radial density distribution of two electrons 
in a cylindrical quantum dot (a = 1). 
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FIG. 12. Singlet (solid line) and triplet (dotted line) density distribution along the z direction 
of two electrons in a cylindrical quantum dot (a = 1). 
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FIG. 13. Energy levels of N e = 1 — 4 electron quantum rings. The parameters of the potential 
are Vq = 10, r\ = 0.5 and T2 = 1. Atomic units are used. 
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2.0 




(atomic units) 



FIG. 14. Density distribution of N e = 1 (solid line), N e = 2 (dotted line), iV e = 3 (dashed 
line) and iV e = 4 (long dashed line) electron quantum rings. The parameters of the potential are 
Vq = 10, r\ = 0.5 and T2 = 1. Atomic units are used. 
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B=0 

B=1.85 




0.0 1.0 2.0 3.0 
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FIG. 15. Density distribution in a single electron quantum ring as a function of the magnetic 
field. The strength of the magnetic field are indicated next to the corresponding density distri- 
bution. The parameters of the potential are Vq = 10, r\ = 2 and r2 = 3. Atomic units are 
used. 
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TABLES 



TABLE I. Comparision of the energies of harmonically confined 2D three-electron system 



(w = 0.2841, hto = 3.37 meV). 
in atomic units. 


The ener 


gies are given 


in meV. Values in parenthesis 


are given 


(M,S) 


SVM 


DIAG | 


] QMC |1 


] QMC || 


QMC 


@ 


(1,1/2) 
(2,1/2) 
(3,3/2) 


26.7827 (2.2582) 
28.2443 (2.3814) 
30.0101 (2.5304) 


26.82 
28.27 
30.02 


26.77 
28.30 
30.04 


26.8214±0.0036 


26.88 
28.35 
30.03 




TABLE 

(, , — o 90 


II. Comparision of the energies of harmonically confined 2D 
Tito = 3.32 meV). rj is the virial factor. 


electron 


systems 


N e 


(M,S) 


QMC 


SVM 




V 


2 


(0,0) 


1.02162(7) 


1.02164 




0.999995 


3 


(IT/2) 


2.2339(3) 


2.2320 




0.999988 


4 
4 
4 


(0,1) 

(2,0) 

(0,0) 


3.7157(4) 
3.7545(1) 
3.7135(6) 


3.7130 
3.7525 
3.7783 




0.999971 
0.999982 
0.999992 


5 


(1,1/2) 


5.5336(3) 


5.5310 




0.999481 


6 


(0,0) 


7.5996(8) 


7.6020 




0.998912 
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TABLE III. Comparision of the energies of harmonically confined 2D three- 
magnetic field(a; = 0.2841). The energies are in meV except for the values in 
are in atomic units, rj is the virial factor. 



electron system in 
parenthesis which 



M,S) 



B(T) 



SVM 



V 



QMC [|11 1 



DIAG § 



(1,1/2) 


0.0 


26.78 


(2.2582) 


0.999991 


26.77 


26.82 


(1,1/2) 


1.0 


26.61 


(2.2442) 


0.999989 


26.60 


26.65 


(1,1/2) 


2.0 


27.69 


(2.3353) 


1.000034 


27.68 


27.74 


(1,1/2) 


3.0 


29.71 


(2.5055) 


0.999987 


29.69 


29.77 


(1,1/2) 


4.0 


32.36 


(2.7283) 


1.000026 


32.32 


32.43 


(1,1/2) 


5.0 


35.39 


(2.9842) 


0.999985 


35.33 


35.48 


(2,1/2) 


0.0 


28.24 


(2.3814) 


0.999992 


28.30 


28.27 


(2,1/2) 


1.0 


27.28 


(2.2998) 


0.999925 


27.33 


27.29 


(2,1/2) 


2.0 


27.67 


(2.3338) 


0.999905 


27.72 


27.69 


(2,1/2) 


3.0 


29.09 


(2.4531) 


0.999954 


29.14 


29.13 


(2,1/2) 


4.0 


31.22 


(2.6324) 


0.999976 


31.26 


31.26 


(2,1/2) 


5.0 


33.79 


(2.8495) 


0.999963 


33.82 


33.85 


(3,3/2) 


0.0 


30.01 


(2.5304) 


0.999999 


30.04 


30.02 


(3,3/2) 


1.0 


28.24 


(2.3817) 


1.000006 


28.27 


28.25 


(3,3/2) 


2.0 


27.97 


(2.3585) 


0.999997 


28.00 


27.98 


(3,3/2) 


3.0 


28.83 


(2.4315) 


0.999999 


28.86 


28.85 


(3,3/2) 


4.0 


30.48 


(2.5703) 


0.999997 


30.51 


30.50 


(3,3/2) 


5.0 


32.63 


(2.7519) 


0.999998 


32.67 


32.66 



TABLE IV. Energies of harmonically confined two-electron system in 3D. 



(L, S, it) 


oj=0M 


w=0.5 


u=10 


(o,o,+) 


0.07921 


2.0000 


32.449 


(1,1,-) 


0.08198 


2.3597 


41.665 


(2,0,+) 


0.08681 


2.7936 


51.338 


(o,o,+) 


0.09696 


2.9401 


52.072 


(3,1,-) 


0.09302 


3.2538 


61.149 


(1,1,-) 


0.10005 


3.3286 


61.504 



37 



TABLE V. Energies of harmonically confined three-electron system in 3D. The SVM is the 
stochastic variational calculation, SM is the shell model and SM-eff is the shell model with effective 
interaction approach (see sec. II D). Atomic units are used. 





SVM SM-eff 


SM 


SVM SM-eff 


SM 


SVM SM-eff 


SM 


(L,S,7T) 


w=0.01 




u=0.5 




co=10 





(1,1/2,- 
(1,3/2,+ 

(2,1/2,+ 
(0,1/2,+ 
(2,1/2,- 
(1,3/2,- 
(1,1/2,- 
(1,1/2,- 



0.181936 
0.182973 
0.184585 
0.191567 
0.198201 
0.193764 
0.193351 
0.199667 



0.181936 
0.182973 
0.184584 
0.191568 
0.187935 
0.193764 
0.193325 
0.199656 



0.181936 
0.182973 
0.184584 
0.191568 
0.187935 
0.193764 
0.193325 
0.199656 



4.013240 
4.310690 
4.366473 
4.467439 
4.717817 
4.794580 
4.805341 
4.960409 



4.013224 
4.310690 
4.366385 
4.467459 
4.717817 
4.794582 
4.797973 
4.957257 



4.013511 
4.310712 
4.366537 
4.467878 
4.717828 
4.794614 
4.798186 
4.957683 



61.138525 
69.972571 
70.315335 
70.853077 
79.490651 
79.860576 
79.890842 
80.793524 



61.138549 
69.972571 
70.315387 
70.853154 
79.490655 
79.860582 
79.890818 
80.793567 



61.139485 
69.972624 
70.315871 
70.854399 
79.490680 
79.860650 
79.891346 
80.794750 



TABLE VI. Energies of harmonically confined four-electron system in 3D. See the caption of 



Table V. 





SVM 


SM-eff 




SM 


SVM 


SM-eff 


SM 


SVM 


SM-eff 


SM 


(L, S, 7r) 




w=0.01 








w=0.5 






w=10 




(1,1,+) 


0.3159 


0.3141 





3141 


6.3492 


6.3490 


6.3502 


91.4466 


91.4459 


91.4496 


(2,0,+) 


0.3177 


0.3188 





3189 


6.3865 


6.3865 


6.3896 


91.6750 


91.6758 


91.6847 


(0,0,+) 


0.3210 


0.3185 





3185 


6.4462 


6.4456 


6.4474 


92.0260 


92.0239 


92.0297 


(0,2,-) 


0.3138 


0.3151 





3151 


6.5875 


6.5875 


6.5879 


99.9068 


99.9041 


99.9053 


(2,0,-) 


0.3198 


0.3181 





3181 


6.7002 


6.6961 


6.6980 


100.5877 


100.5875 


100.5930 


(1,1,-) 


0.3240 


0.3195 





3195 


6.7196 


6.7093 


6.7105 


100.6478 


100.6199 


100.6235 


(i,o,-) 


0.3278 


0.3251 





3251 


6.7961 


6.7935 


6.7963 


101.0946 


101.0740 


101.0813 


(1,1,-) 


0.3408 


0.3232 





3232 


6.8448 


6.8153 


6.8169 


101.3253 


101.2220 


101.2270 


(2,2,+) 


0.3223 


0.3212 





3212 


7.0385 


7.0202 


7.0205 


109.5179 


109.5156 


109.5162 


(1,2,+) 


0.3264 


0.3313 





3313 


7.0702 


7.0706 


7.0719 


109.7618 


109.7612 


109.7638 
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TABLE VII. Energies of harmonically confined five-electron system in 3D. See the caption of 



Table V. 





SVM 


SM-eff 




ni\ it 

SM 


SVM 


SM-eff 


SM 


SVM 


rn\ it nr* 

SM-eff 


SM 


(L, S, tt) 




cj=0.01 








u=0.5 






L0 = W 






(0,3/2,-) 


0.4804 


0.5141 





5165 


8.9963 


8.9979 


9.0032 


123.357 


123.3539 


123 


3682 


(2,1/2,-) 


0.4858 


0.5175 





5203 


9.0567 


9.0526 


9.0588 


123.749 


123.6960 


123 


7129 


(MA-) 


0.4880 


0.5186 





5211 


9.0954 


9.0919 


9.0988 


123.949 


123.9287 


123 


9482 


(1,3/2,+) 


0.4869 


0.5318 





5359 


9.3110 


9.2969 


9.3024 


132.320 


132.1385 


132 


1523 


(0,1/2,+) 


0.4931 


0.5450 





5525 


9.4443 


9.4355 


9.4458 


133.045 


132.9021 


132 


9265 


(1,3/2,+) 


0.5108 


0.5472 





5537 


9.7104 


9.4701 


9.3692 


133.223 


133.1205 


133 


1427 


(2,1/2,+) 


0.4950 


0.5357 





5406 


9.3582 


9.3528 


9.3599 


133.471 


132.487 


132 


5026 


(0,1/2,+) 


0.5267 


0.5561 





5644 


9.8766 


9.5866 


9.5990 


134.204 


133.8337 


133 


8658 


(2,5/2,-) 


0.4829 


0.5232 





5253 


9.5919 


9.5891 


9.5914 


140.973 


140.9054 


140 


9105 


(0,5/2,-) 


0.4882 


0.5306 





5336 


9.6626 


9.6618 


9.6648 


141.270 


141.2692 


141 


2762 



TABLE VIII. Energies of harmonically confined six-electron system in 3D. See the caption of 



Table V. 





SVM 


SM-eff 


SM 


SVM 


SM-eff 


SM 


SVM 


SM-eff 


SM 


(L,S,7T) 




o;=0.01 






^=0.5 






u=10 




(1,1,+) 


0.703 


0.797 


0.815 


12.038 


12.064 


12.079 


157.701 


157.415 


157.451 


(2,0,+) 


0.743 


0.801 


0.819 


12.080 


12.101 


12.118 


157.910 


157.643 


157.681 


(o,o,+) 


0.714 


0.805 


0.822 


12.128 


12.159 


12.178 


158.080 


157.991 


158.034 



TABLE IX. Energies of harmonically confined seven-electron system in 3D. The SM is the shell 
model calculation. 





SM 


SM 


SM 


(L,S,7T) 


w=0.01 


u;=0.5 


w=10 


(1,1/2,-) 


1.063 


15.390 


193.055 


(1,3/2,-) 


1.084 


15.934 


209.998 


(3,5/2,-) 


1.087 


15.960 


210.134 


(0,1/2,+) 


1.104 


15.672 


201.377 


(2,3/2,+) 


1.108 


15.707 


201.575 


(3,3/2,+) 


1.113 


15.749 


201.802 
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TABLE X. Energies of harmonically confined eight-electron system in 3D. See the caption of 



Table IX. 







bM 




niv it 

bM 


bM 


(L, S, 7r) 




w=0.01 




UJ=0.& 


LO=10 


(0 o +) 




1.412 




19.038 


230.219 


(0,2,+) 




1.448 




19.650 


247.204 


(2,1,+) 




1.448 




19.653 


247.212 


(1 1 -1 




1.475 




19.430 


238 771 


(3,1,-) 




1.479 




19.456 


238.915 


(2,0,-) 




1.483 




19.491 


239.131 


(2,1,-) 




1.484 




19.496 


239.139 




TABLE XI. Properties 


of harmonically confined 2D systems. 




N e (M,S) 






u = 0.01 


uj = 0.5 


uj = 10 




(H) 




0.0738 


1.659 


23.652 




(T) 




0.0092 


0.443 


9.29/ 


2 (0,0) 


(Vcoul) 




0.0369 


0.516 


6.6 ( 2 




(V^con) 




0.0277 


0.701 


lU.yoo 




V 




0.9999998 


0.9999995 


0.9999998 




(H) 




0.176 


3.573 


48.365 




(T) 




0.016 


0.822 


lO.ZOD 


3 (1,1/2) 


(Vcoul) 




0.096 


1.286 






(V^con) 




0.064 


1.465 


oo ion 




V 




0.9999972 


0.9999984 


1 1 1 1 1 1 1 1 1 1 1 1 . i 

0.9999981 




(H) 




0.317 


5.863 


74.979 




(T) 




0.018 


1.137 


26. 836 


4 (0,1) 


(Vcoul) 




0.186 


2.391 


I yj ICO 

14.1o3 




(V^con) 




0.112 


2.335 


33.981 




V 




0.999812 


0.999921 


0.999942 




(H) 




0.515 


8.670 


104.642 




(T) 




0.0196 


1.421 


34.931 


5 (1,1/2) 


(Vcoul) 




0.339 


3.874 


23.168 




(V^con) 




0.159 


3.376 


46.543 




V 




0.9992 


0.9995 


0.9991 
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TABLE XII. Properties of harmonically confined 3D systems. 


N e (L,S,7r) 




uj = 0.01 


u = 0.5 


u = 10 




(H) 


0.0792 


2.0000 


32.4486 




(T) 


0.0121 


0.6644 


14.4412 


2 (0,0,+) 


(Vcoul) 


0.0366 


0.4474 


2.3776 




(V con ) 


0.0304 


0.8881 


15.6299 




V 


0.999999 


0.999999 


0.999999 




(H) 


0.1819 


4.0132 


61.1385 




(T) 


0.0192 


1.1507 


26.0867 


3 (1,1/2 -) 


(Vboul) 


0.0957 


1.1411 


5.9763 






0.0671 


1.7214 


29.0755 




V 


0.999991 


0.999995 


0.999999 




(H) 


0.3161 


6.3502 


91.446 




(T) 


0.0229 


1.5853 


37.371 


4 (1,1,+) 


(Vboul) 


0.1770 


2.1174 


11.132 




(Von) 


0.1163 


2.6475 


42.943 




V 


0.999821 


r\ a Ann a 1 

0.999891 


r\ r\r\r\r\-t r\ 

0.999912 




(H) 


0.48041 


8.9963 


123.36 




(T) 


0.02501 


1.9786 


48.283 


5 (0,3/2,-) 


(Vcoul) 


0.27881 


3.3562 


17.808 




(Vcon) 


0.17660 


3.6615 


57.266 




r) 


0.99812 


0.999671 


0.999781 
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